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NOMENCLATURE 



Cp Specific heat of steel 

D Diameter of ingot 

E Heat of fusion of steel 

h i Surface coefficient of heat transfer, ingot to mold 

farj Surface coefficient of heat transfer, Mater side 

1 Mold wall thickness 

L Mold length 

q Heat transferred, Btu/br.,ft2 

Q Heat transferred, Btu/hr. 

ri radius of ingot 

V2 radius of mold, water side 

th average surface temperature of ingot within mold 
i>2 average water temperature in mold 

mean temperature of ingot cross-section at mold exit 
T p pouring temperature 

U Composite heat transfer coefficient 

velocity of ingot 
thermal diffueivity 



density of steel 
hwC>/^ Susselt number, water 
De »/K Reynold^ number, water 
Prandtl number, water 
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SUMHAEY 



Analyst a of continuous casting involves many phases. Wo of which 
are described in this thesis* The first investigation consists of 
determining the temperature distribution in a moving Ingot of metal 
which is solidified and cooled sufficiently in a water cooled mold 
to permit the continuous processing of the resulting billet* In 
the second phase there has been established a relationship for maximum 
water flow rates for heat transfer and their association with cor- 
responding coefficients of heat transfer* 

The analysis has shown that the rate of heat transfer is primarily 
a function of the surface coefficient of heat transfer between the 
Ingot and the mold* To increase this rate some means must be provided 
to increase this coefficient. Addition of oil along the boundary has 
been suggested which would not only increase this surface coefficient 
but will also decrease the frictional drag* 

The results, although only approximate, are provided as a guide 
for future work. 
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INTRODUCTION 



Generally, the analytic mode of approach to the solution of a 
problem also requires experimentation. However, lack of facilities and 
time limitations precluded experimentation for verification of any 
results in this thesis* 

The primary endeavor in this investigation is to assist the 
experimenter in his work by forming limits of practicability and 
providing mathmatical solutions to important phases of the problem. 
Results are only approximations, consistent with assumptions necessary 
to perform the operations, but they can be used as a guide for future 
work. 

Treatment of the entire problem of continuous casting is very 
long and complex, and to the best knowledge of the writer no work has 
been done analytically along these lines, either in full or upon any 
phase. No information is available for many of the problems; however, 
with work being done at the present time on properties of metals at 
temperatures at and near the melting point, some hope is held that in 
the future information will be available concerning these properties. 

The process involves pouring molten steel (with a minimum amount 
of superheat to allow good fluidity) into a vertical, water cooled, 
metal mold. At the start of the process the mold is plugged at the 
lower end until sufficient freezing takes place to allow adequate 
strength to permit pulling the plug and continuously processing the 
resulting billet* 

An analysis involves investigation of various phases which are 
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outlined below and discussed later In this paper* 

(a) Temperature distribution at any point in the ingot. 

(b) Investigation of limiting water flow rates, and 

(c) Prediction of accompanying coefficients of heat transfer 
between the ingot and the mold. 

(d) Possibility of frictional drag exerted on the mold by 
the ingot. 

(e) Maximum velocity of ingot travel. 

(f) Load capacity of the ingot. 

(g) Determination of ingot breaking point. 

(h) Investigation of possible surface boiling on water side 
of Idle mold. 

(i) Investigation of surface tension and surface film effects. 
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HISTORY 



The greatest amount of work on continuous casting of steel has 
been done by Edward R. Williams and the Williams Engineering Company, 
His solutions to the problems hare been of an experimental nature, 
and to the best knowledge of the writer, no attempt has been made 
toward an analytical investigation* 

The Williams process accomplishes the necessary heat removal by 
providing a very thin mold wall of high conductivity material (com- 
mercial drawn copper or brass tubing) through which heat can rapidly 
pass from the ingot, and by removing this heat extremely rapidly from 
the outside of the mold wall by means of water or other cooling fluid 
passed at high velocity along the surface* By removing the heat 
rapidly from the outer surface of the mold vail, a high temperature 
gradient is produced in the mold wall, which permits fast removal of 
heat from the ingot* It is claimed that this instantly freezes the 
molten metal into a skin which shrinks away from frictional resistance 
and is sufficiently strong to prevent fracturing as it is moved along 
the mold wall. 

The cooling water in the mold passes through passages which have 
a very small clearance, the velocity being quite high. It has been 
found that a relatively small amount of water is sufficient to satis- 
factorily remove the heat from a steel ingot to enable it to be poured 
and withdrawn at a speed of seven or more feet per minute* 

As the ingot emerges from the bottom it passes through a series 
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of sizing rolls that prevent bulging of the walls by ferrostatic pres- 
sure from the molten interior, "Water sprays surround these sizing 
rolls and the ingot is quite cool when it gets down to the pinch rolls, 
Williams frequently uses mold lubrication to increase speed of cas- 
ting, The lubricant is introduced at the top along the mold wall, 
the quantity being controlled by the speed of the pinch rolls drawing 
the continuous billet out. 

The Williams machine has been given no steady trial over a long 
run on steel; however, the experimental ingots produced have, for the 
most part, satisfactory structure, 

A number of makers of carbon steel who habitually handle small 
heats have shown interest in the process. It is also conceivable 
that those makers of high alloy steels and certain types of stainless 
steels may see in this process a possible means for eliminating present 
high scrap losses and a more flexible method of handling small tonnage 
orders. 
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DETERMINATION OF THE TEMPERATURE AT ANT POINT IN THE INGOT 



la the following Investigation certain assumptions are made to 
simplify the problem and bring it within the realm of possible solution. 
These assumptions are as follows! 



(a) Thermal conductivity does not vary with tem- 
perature. 

(b) No recognition of the heat of fusion of steel 
Is made in the temperature distribution in 
the ingot. 

(c) Convection in the molten metal is neglected. 



It is realized that the first two assumptions will cause appreciable 
error in the final results; however, unless these assumptions are made 
the problem would be beyond the scope of this paper. The assumption 
of the absence of convection in the molten steel is justified since it 
will become "mushy* upon entrance into the mold and resist any ap- 
preciable turbulence. 

For purposes of this analysis a cylindrical ingot is used with 
all results in cylindrical coordinates. Consider first an elemental 
volume as follows* 
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The partial differential equation expressing the heat conduction 



for this problem lss (refer to Appendix n A* for derivation) . 
b'f’ 2 " 



6? 1 



bt 



~ O 



vhero c*f la the •thermal dlffasivity* of the material. 



~r -- 4^9 



The thermal diffusivity la a property of the substance in -which 
heat is transferred and stored. It Is the ratio of a quantity which 



is proportional to tho heat conduction (numerator) and to a quantity 
which is proportional to the heat storage (denominator). Heat units 
cancel and the dimensions become ft^/hr. 

Assuming it is possible to express the solution of the above 
equation as the product of 3 functions, one of which is a function of 



■r* alone, one a function of 'tfl* alone, and the other a function of 
•a* alone, then 



t- RV, £(<!>) 2 fe> 



Tho solution (as shown Appendix «B") is 















where ■a® is a constant 

Note* Since the temperature distribution is a function of 
B r* and *a" only, terras including *<$)» drop out. 
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The Ingot is considered to be a semi-infinite cylinder with the 
following boundary conditions* 

1* t- = 0 , o<r<.R 

2* t^-o , % r 

3* f ~R«s> , o<«<-<» , r = ^ 

From the second boundary condition it is seen that for 
at ? = , C 3 must be zero* The equation then reduces to 

t = 

Solution for the above is made by assuming that the temperature n t n is 
the sum of two solutions, that is 



Solutions for t. and t t are taken as modified forms of equations 
31 and 32, pp 192, Carslaw and Jaeger, Conduction of Heat in Solids or 



and the a ^ are toe positive roots of JjRa^ s 0, 

The total solution then becomes the sum of equations (1) and (2) above* 
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Reducing the latter part of equation (2) to more usable form, we have. 



This farther reduces to 



2 * |p(p) c K " ( f’tj(s - ,,( '^' 1 



The temperature distribution in the ingot thus becomes 

-fc 

^(ra,*) r -V 



w n is 

2 V 3 0 (jraJ / 0 

t- H. m.> e i 



Z 



£(c<0 



Z j F(f) -ymK K n fz-^)ci|£ 



R ^ 3<fca J U 

0-1 ' •— 

A(f) e”^ - /f-^le^'dp 



For actual temperatures the functions f(r) and f must be 
supplied from physical aspects of the problem. Consider, for steady 
flow, the following assumptions j 

f(r) z Tp (pouring temperature) 

F( ) » Tpe^where L is the length of the mold. 
Substituting these values and integrating we have 
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further simplification, 


, this reduces to 


Up V ^a*) * 
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Since K n is negative for all values of this vill converge very 
rapidly. 
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INVESTIGATION OF MAXIMUM COOLING WATER FLOW RATES, RELATION TO HEAT 
TRANSFER COEFFICIENT, AND EXPRESSION FOR MOLD LENGTH 



For this phase of the problem consider the magnitude of the 
heat transferred from the ingot to the cooling water with the as- 
sumption that the velocity of the ingot is such as to cause final 
freezing at the instant the ingot emerges from the mold. Since the 
mass rate of flow of the metal is constant, the heat transferred in 
unit time for a mold of length L is equal to 



Q a AChin - h out ) 


Btu/sec 




. *TT t> 1 / 

where m — - ^ 


ii it ii 


velocity, ft/hr 
diameter ingot 
density steel 


“ hin . H + C f Ip 


Ip s 


pouring temperature 


H = 


heat of fusion 




Cp r 


Specific heat 




T - 


mean temperature of 
ingot cross-section 
at mold exit 



Substituting these values into the above equation, the total heat 
transferred is 

Q - [4 +C ( ,C V t')>./u 

Dividing through by tIDL, we have 
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Since p q" and "V^" are Independent there is obtained an expression 
for the length of the mold in terms of the heat transferred. 
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Prom the basic law of heat conduction for a hollow cylinder. 



the heat transferred from the ingot to the cooling water is 

- U. L 2 TT<~ # 
dr 



and by integration within the limit* r^ and r 2 
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J2.vi r - {?* r t 



where 

ri s radius of ingot 
r 2 = radius of mold, water side 
tl • average surface temperature of ingot within 
mold 

to Z average water temperature in the mold 



The term "U" is a composite or "overall" coefficient of heat transfer 
and is defined as follows I 

_i_ j_ , 4 4- 

u" ^ + ^ 

where 

hy s surface coefficient of heat transfer, water 
side 

1% z thermal conductivity of mold 
h^ s surface coefficient of heat transfer, ingot 
to mold 

1 z mold wall thickness 



from this, it is seen that 
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Substituting this into the equation of heat conduction and solving 



for h; , we have 



- "tiw (i* r % - Q 

^ ^ ('tw+Hw)(fvir | -io f^-r (v^x) 



For liquids with Beynold*s numbers exceeding 10,000 the expression 
for h w is given tgr the Colburn equation 

(equation Ud, pp 168, McAdams, Heat Transmission ) 

(^-)= Nms^e\\ Ma.mVer 

~ Revjvio\i* NlvxYTkWec 

I^J^j - Vra^dti fOuwWc 



Applying this data to a specific problem, consider the fol- 



lowing datat 



Diameter of ingot 
Pouring temperature, T 
Length of mold, L 
Freezing temperature 

7 

Heat of fusion of steel 
Specific heat steel 
Mold thickness 
*wA 



h inches 
2800°F 
3 ft. 

2760°F 
2560OF 
65 Btu/# 

0.1075 Btu/# 
0.5 inch 
1650 B/hr.ft.F 
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Frost this 



Q = 11.1 x 10j> Btu/hr 
q : 3.51» I 1( P Btu/hr,ft 2 

Using this value of "q" in the equation for the determination of 
the surface coefficient of heat transfer (ingot to mold), h-^ can be 
determined. With this value and values for h w and k w /l, an "overall” 
coefficient for heat transfer for varying Reynold’s numbers is 
calculated. Curves of h^, h w , ky/l and U versus Reynolds numbers 
are shown plate I, pp iM. 

The results of these curves show that the overall coef- 
ficient of heat transfer is dependent upon the surface coefficient 
between the ingot and the mold, and for this particular set of con- 
ditions, water flow rates in excess of Reynolds number 5 x 10^ would 
be impractical. Such determinations can be made for varying sets 
of conditions giving a family of similar curves, and from these flow 
rates, mold length, and surface coefficients of heat transfer can 
be predicted. These results are not conclusive as much depends 
upon reasonably accurate determination of temperatures within the 
ingot. It does, however, show that maximum cooling water flow rates 
are dependent upon the amount of heat that can be transferred from 
the ingot to the mold, which in turn is dependent upon the correspond- 
ing surface coefficient. 
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Coefficient of Heat Transfer 




CONCISIONS 



1. As pointed out previously the solution for the temperature 
distribution is only approximate because the heat of fusion of the 
metal has been neglected in the derivation. The solution, however, 
is useful in heat transfer problems for determining temperatures in 
cylindrical bodies moving with constant velocity, and also in 
predicting temperature distributions in continuous casting problems. 
Solutions for different velocities can be modified by proper selection 
of the functions f (r) and which can be fairly accurately pre- 
dicted from experimental data* 

2. In the investigation of maximum water flow rates, and their 
relation to surface coefficients of heat transfer between ingot and 
mold, it is shown that the over-all coefficient of heat transfer is, 
in the limit, a function of the surfaoe coefficient between the billet 
and the mold. From the basic law of heat conduction and for a given 
3lze ingot, it is concluded that to increase the heat transferred some 
moans must be provided to increase this surface coefficient. In the 
Williams process mold lubrication is used, which, in addition to re- 
ducing friction a resistance, increases the surface coefficient of 
heat transfer by making the heat transfer across the space more nearly 
that of pure conduction. It Is recommended that for the lubricant 

an oil of high thermal conductivity be used - the optimum between 
the lubricating effect and the increased heat transfer coefficient is 
left for future investigation. 
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After determination of the maximum heat transfer rate, the 
optimum length of the mold can be determined from the relation 



It is to be noted here, however, that the above equation is based on 
the assumption that final freezing takes place just as the ingot 
emerges from the mold* More realistically, this final freezing point 
takes place some distance below the mold* That is, the mold length 
need be only sufficient to establish a "skin” around the molten metal 
which is strong enough to prevent rupture either from excessive load 
or from ferrostatic pressure from the molten interior. Investigation 
of minimum mold length is recommended for further work. 
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RECOMMENDATIONS FOR FURTHER WORK 



Since the technique of continuous casting is complex, and 
since time limitations precluded analysis of all phase?, the following 
recommendations are made for considerations 

1. Frictional drag exerted on the mold by the ingot 

Although it is claimed that the molten metal instantaneously 
freezes upon entry into the mold and shrinks away from frictional 
resistance (The Iron Age, April It Mid April 11, 19l;0), it seems highly 
improbable. One of the laws of sliding friction is that the friction 
is proportional to the force pressing the surfaces together. This 
force is increased by the ferrostatic pressure from the molten 
interior, which, it is believed, will overcome the shrinking due 
to freezing of the metal. Adding a lubricant to decrease the coef- 
ficient of friction will increase the maximum load capacity and hence 
the maximum allowable velocity. 

Investigation of the mechanical friction is simplified be- 
cause, within limits, sliding friction is independent of the area of 
contact, and (except at start) of the speed of relative motion. 

The problem, therefore, resolves itself to a balance of forces be- 
tween the ferrostatic pressure and the shrinking due to cooling. 

2. Maximum Velocity of ingot travel. 

The velocity of the ingot varies directly as the amount of 
heat which can be satisfactorily removed from the steel. With suf- 
ficient heat removal from the mold a skin will be developed such that 
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the Ingot can withstand the load imposed by the pinch rolls* The 
maximum velocity will be determined by the maximum heat removal 
consistent with strength characteristics of the skin of the billet* 

3* Load capacity of ingot and determination of 
ingot breaking point 

The load capacity of the ingot is a function of the load 
capacity of the skin formed in the mold* The load will be caused 
by surface tension and surface film effects of the steel and friction- 
al resistance to the passage of the ingot. This load capacity will 
reach its maximum when such velocity is reached that the maximum 
heat removal will be below that necessary to develop sufficient skin 
thickness and strength* Since the ferrostatic pressure increases 
down through the mold, the frictional resistance will also increase, 
and hence some point will be reached where the load exceeds the 
ultimate strength of the billet. Experimental evidence has shown 
that the breaking point will occur some distance down the mold, but 
the exact location has not been determined. 

i* Investigation of possible surface boiling on the water 
side of the mold 

With very high rates of heat transfer and high temperature 
gradients produced in the mold wall there are excellent possibilities 
of surface boiling in the cooling water. With high feed rates, a 
small fraction is vaporized, and the boiling-side coefficients may- 
be roughly twice those predicted from straight water cooling. These 
increased coefficients may be due to increased turbulence in the film 
near the wall. Investigation along these lines can be very instru- 
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mental in increasing heat transfer rates* 

•>* Surface Tension and Surface Film Effects 
A high value of surface tension has the effect of increasing 
the pressure required to cause the metal to enter and flow along a 
narrow channel} but the direct effect of the true surface tension 
is almost overshadowed by the effect of surface films , the influence 
of which is equivalent to that of a greatly increased surface tension* 
It has been found that surface oxide films affect the flowing power 
of metals with the apparent surface tension of some metals increased 
three times* Steel has a very high value of Burface tension* 

(Briggs, Metallurgy of Steel Castings) 

Both effects will add to the resistance to flow and decrease 
the maximum load capacity of the ingot. Although it is thought to 
be considerably small compared to frictional resistance, further 
investigation should be initiated to determine its effects* 
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APPENDIX A 



Equation of Conduction for a Moving Cylinder 




Beat conducted through "r" face is 

rJ<pdzJ>T - - - - - dQ xr 

Heat conducted through the *r -»■ dr* face is 

kdcfdz—rdr + kd<P</z-^JrJr -t-Udd* r^drdt - - - dQ 
br dr * r 

Heat conducted through *<j? w face 

k|«' lrd ' jT d °'« 

Note* Temperature gradient in cp direction is the 
change of temperature with respect to the 
distance represented hy the angled which 
distance also depends on the radius r. 

Heat conducted through * cP + dp * face 

J E fcfap) ^ ^d?) Cfr<f) rc Kl^ r 
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Heat conducted through "a" face 



t, ^ r d4>drclr 

N it 



<*?« 



rdcpdrdnr 4- kd<Pd<- r ^7 



Heat conducted through, "z +dz B face 

hi 

hz 

or |< r<J<pdrdv 4- ko/^rdir^dr d#.,,* 

Heat gained by the elemental volume due to conduction is equal to 

r^<v I d<3 ‘ 

which equals 

f y^ddlirdih + dT + A^^dkJrd.tolT V f^r 



k 



Heat stored in the differential body 

<-}Q 5 = ^ C f (cdcjdzdr^cl'f 

However the temperature increment of a particle of the fluid in 
time interval dT becomes . , 
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dividing through b* dT 

di \, 






2 62 



Jyt_ 

6T 



iiace we are consideriit velocity in the *-41rection end a condition 
of steady state, this reduces to 

From this, then, the Heat stored In the body become 

Aq 3 ^ ^ c r (^cPdrcJe) v t -| | dr 

Since no heat energy ie being developed in the body, the law of the 
conservation of energy reduces to 

Meet stored s heat in - heat out 

dQ 3 - <i£) x - d Q 2 

Substituting and dividing through by rd<P^l t'diclT 

4-' 



Q r j Vl-vT*^ > Vi . Vf| 
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introducing the terra "thermal difXtoaivlty 1, 



^ = k / 



'( ,c v 



(The thermal diffusivlty is a property of the substance in which 
heat ie transferred and stored)* 

Substituting tola value, the equation becomes 

Xi . J- Vi Vi 



^ * v * V- ^ 



_ ^ = o 

iz 1 °< 
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APPENDIX B 



SOLUTION FOR THE EQUATION OF CONDUCTION - TEMPERATURE DISTRIBUTION IN 
INGOT 



A seme it is possible to express the solution of the dif- 
ferential equation (Appendix A) as the product of three functions, 
one of which is a function of "r" alone, one a function of 
alone, and the other a function of ”z" alone, then 

t = R( r) S>(<$) %,(£) 

The equation takes the fora, 

Z 4 -^ 2 , $> R' <- 7 , ^'R + Z"& R-^ - o 

dividing through ty R<£>2« * transposing and equating to a constant, 

ve have 

R + r r + r 1 & Z z 

From the left side of the equation, write 

r z — +■ +• ~ 

R R ® 

which gives 



f 



• a o) 
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the solution of which is 



<4> - +■ C z Cos p (ft 
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Equation (2) also gives 



r*R ’ +■ r R' 4-(rV-ip x ) R = o (5") 

let x : ar 

then since . 

dR dR dx an d a 
Tr dt dr dr 

equation (5) becomes 

mW R= W = RqO 

This equation is Bessel's differential equation (Hildebrands Advanced 
Calculus of Engineers, pp 165-167), and the general solution for "p* 
an integer can be written 

£>= C, 3 r (x) +C z K p (*) 0) 

where Kp (x) is given as follows (Woods, Advanced Calculus, pp 288) 

lsf>W = 5*> W + fix) 

. ft, , r n 

u) W Ho = 2. l_ 2-r^" n i 2 . L z r 1 "nW! * J 

h-o * h-o ' 

Since equation (7) is applied to a solid cylinder, Bessel's function 

Kp(x) becomes infinite at x : 0 (in 0 s -»o) , and it follows that 

C 2 = 0. Thus for the solid cylinder, the general solution is taken 

as: 

£>= C = C Jp(,ar) (.8) 
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( 3 ) 



r .. 



©) 



Taking the remaining equation in our solution 



Using the operator method of solution, the equation becomes 



The solution, therefore, of equation (9) becomes 






= C>e 



V c 



.<2 




(10) 



O') 



Hence for the cylinder the solution of the differential equation is 
given ty substituting equations (J|) , (8), and (ll) into equation 



t= R(y) 

iWn t= C3^hr)|c, Svwvp(l) vC-jCop «p] fce’+c/ 1 ] 



Since the temperature distribution is a function of "r B and "z" only, 
then, p * 0, and the solution of the equation becomes 



t= 
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